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Abstract. The large number of low- lying states of d- and /-shells supports a variety of 
order parameters. The effective dimensionality of the local Hilbert space depends on the 
strength, and kind, of intersite interactions. This gives rise to complicated phase dia- 
grams, and an enhanced role of frustration and fluctuation effects. The general principles 
are illustrated on the example of the effect of a magnetic field on quadrupolar phase 
transitions in some Pr-based skutterudite compounds. 



1. Introduction 

Transition metal and rare earth compounds show a rich variety of collective 
behavior: various kinds of ordered phases as well as strongly fluctuating 
states (spin and orbital liquids). The basic reason is that d- and /-shells 
have a relatively large number of low-energy states. Crystal field splitting 
usually reduces this number (the dimensionality D of the local Hilbert 
space) considerably below the free ion value, but complicated physics can 
arise even from T> = 3 or 4. 

Let us briefly consider some examples. T> = 3 is, in one interpretation, 
the case of S = 1 spin models which turned out to have unanticipated 
phases like the spin nematics [1]. A different realization is offered by in- 
teracting /-electron models based on the low-lying quasi-triplet of Pr ions 
in PrBa2Cu307_<5 where the nature of Pr ordering is still an open issue 
[2, 3]. A literal realization of T> = 4 is offered by the T% ground state 
of Ce ions in CeBg which has a rich phase diagram [4]. Alternatively, we 
may think of the fourfold quasidegeneracy arising from the combination 
of twofold spin degeneracy with twofold orbital degeneracy, which is the 
simplest model of cf-electrons which is capable of supporting either spin 
or orbital order, or a combination of both. Twofold orbital degeneracy 
may occur in cubic, tetragonal, and hexagonal environments, and there 
are quite different versions of the e g -model. In the cubic e g model, orbital 
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momentum is completely quenched, but orbital order may still break time 
reversal invariance in the octupolar phase [5]. In contrast, the trigonal 
e g states sustain permanent orbital momentum along the threefold axis; 
such a model should be relevant for understanding the complex behavior of 
BaVS3 [6, 7]. The possibilities of V > 4 models are largely unexplored, but 
we should mention the V = 6 t2g-models of LaTiC-3. There is an intricate 
relationship between orbital ordering and spin ferromagnetism which would 
be remarkably difficult to explain without due consideration of the orbital 
degrees of freedom [8] . 

As implied above, the definition of T> is not straightforward: it is usually 
higher than the degeneracy of the ground state level (which would often 
be small because of small low-symmetry components of the crystal field), 
but it is not so high as the free-ion value. In CeB@, the fourfold degenerate 
Tg, level is well separated from the higher-lying Tj which still arises from 
the Hund's rule ground state; in PrP^CusC^-^, a low-lying doublet and a 
singlet can be lumped together to give a quasi-triplet which would have Fx, 
character if the symmetry were cubic; but it is, in fact, only tetragonal. In 
any case, the relevant dimensionality T> = 3 is much smaller than 9 which 
would be the Hund's rule value. It depends primarily on the strength of 
intersite interaction, which splittings should be considered small. 

The highest symmetry of the the V = 3 models would be SU(3), etc. 
Clearly, the exact realization of a high-symmetry model is more than im- 
probable, and if it were really required, we should forget about it. However, 
there are indications that the domain of influence of such a seemingly arti- 
ficially high symmetry point in parameter space extends over a substantial 
portion of the phase diagram [7]. It stands to reason that SU(D) models 
(which have symmetries connecting spin and orbital axes in Hilbert space) 
are more quantum fluctuating than the pure SU(2) spin models: there are 
more transverse directions to fluctuate to. For instance, the SU(4) model 
on the triangular lattice has a plaquette resonating ground state (the SU(4) 
version of the resonating valence bond idea [9]), and the influence of this 
spin-orbital liquid state may extend to physically relevant regions of the 
parameter space [10]. 

Limitation of space forbids us to present more than one concrete exam- 
ple of spin-orbital models. We consider Pr-filled skutterudites, in particular 
PrFe4Pi2 which we model as a V = 4 system. 

2. The case of PrFe 4 Pi 2 

Pr-filled skutterudites show varied behavior: PrRu4Pi2 has a metal-insula- 
tor transition [13], PrOs4Sbi2 is thought to be an exotic superconductor 
[14], while PrFe4Pi2 remains a normal metal in the entire temperature 
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range studied so far. Our interest lies in PrFe4Pi2 which in a certain 
parameter range can be characterized as a heavy fermion system with 
exceptionally high electronic specific heat [11]. PrFe4Pi2 has a phase tran- 
sition at Tt r ~ 6.5K to an ordered phase which had first been thought to 
be antiferromagnetic, but mounting evidence indicates that it is, in fact, 
antiferroquadrupolar (AFQ) [11, 15, 12]. 

Our purpose is to model the AFQ transition by a crystal field model 
of the 4/ electrons. We assume that Pr is trivalent (4/ 2 ). This can not 
be literally true, since the driving force of the formation of a heavy band 
is presumably the admixture of other valence states. However, high-field 
studies show that the heavy fermion state competes with AFQ ordering [11], 
so a localized /-shell description should be acceptable within, or adjacent to, 
the AFQ phase. Besides, at low fields H and high temperatures T, thermal 
dehybridization acts to obviate the need to consider interband coherence 
effects. 

2.1. THE CRYSTAL FIELD MODEL 

The J = 4 manifold of Pr 3+ is split by the approximately cubic crsytal 
field into the T± singlet, the T3 doublet, and the T4 and V*, triplets 1 . Group 
theory does not tell us the sequence of the states, but fitting the measure- 
ments narrows the choice. Analyzing the anisotropy of the magnetization 
curves, it was concluded that the likely possibilities are: a T\ ground state 
and a low-lying exited state T4 (the Ti-^ scheme); or the ri— T$ scheme; 
or the T3-F4 scheme [11]. Similar schemes were suggested for PrOs4Sbi2 



The assumption that the low-T ordered state is AFQ, seems to speak in 
favour of the T3— T4 scheme, since then the ionic ground state T3 posesses a 
(permanent) quadrupolar moment. It was also pointed out that the choice of 
the T3 ground state is consistent with a symmetry analysis of the strucrural 
distortion accompanying the AFQ ordering [12]. This latter argument relies 
only on the assumption of the T3 ground state, and does not consider the 
effects of the low-lying excited state. Here we show that the assumption of 
the Ti-r4 scheme is also capable to account for most of the observed static 
properties of PrFe4Pi2- 

We now discuss the consequences of assuming a Ti-^ level scheme. 
Since the singlet ground state 



1 The symmetry group is really not Oh, but the tetrahedral Th- We nevertheless use 
the cubic classification, which is an approximation at zero field, but when H 7^ 0, the 
symmetry will be in any case substantially lowered. 



[14, 16]. 




(1) 
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does not carry any kind of moment, the ordered quadrupolar moment has to 
be induced by intersite interactions, assuming that the local Hilbert space 
contains also the triplet 

|r 4 + ) = J{|3> + |-3> + >/7(|l> + |-l»} (2) 
|r2> = ^(|-4)-|4)) (3) 

|r 4 -) = | {|3> _ | _ 3 > + (I - 1> - |i>)} , (4) 

where we have chosen the basis of quadrupolar eigenstates. Choosing the 
energy of (1) as the zero, the states (2)-(4) lie at the level A. 

The possible moments in the four-dimensional local Hilbert space spanned 
by (l)-(4) are given by the decomposition 

(Ti e r 4 ) (ri © r 4 ) = 2Ti + r 3 + 3r 4 + r 5 . (5) 

Evidently, the system could support either dipolar (T 4 ), or either of two 
kinds of quadrupolar (T3 or T$) order 2 . The quadrupolar order parameters 
are the same that appear in the decomposition of a purely T 4 system 

r 4 ® r 4 = ri + r 3 + r 4 + r 5 , (6) 

i.e., they are not sustained by inter- level matrix elements 3 . - It may be 
of some interest to mention that the Ti-r 4 scheme does not offer the 
possibility of octupolar order (but the T3-r 4 scheme would). 

2.2. THE EFFECT OF EXTERNAL MAGNETIC FIELD 

Our decomposition (5) allows to seek dipolar and/or quadrupolar ordering 
in the system. Experiments give the clue that we should, in fact, look for 
(antiferro) quadrupolar order. We may rather arbitrarily assume that it is 
of the T3 kind 4 , i.e., the possible order parameters are 0\ = 3J 2 — J(J + 1) 
and 0\ = J% — Jy. Furthermore, since the total energy expression for a pair 
of sites has only tetragonal (as opposed to cubic) symmetry, we need not 
assume that the 0% and 0\ couplings would be equal, and we may seek, 

2 This is a classification of order parameters which can be defined purely locally. Q^O 
order needs further discussion. 

3 In contrast, matrix elements between Ti and T4 would bring extra possibilities of 
dipolar ordering. This does not seem to be relevant for PrFe4Pi2- 

4 Assuming O xy -type Fr, ordering would give similar results. 
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say, C|-type order. Using a mean field decoupling, our problem would be 
rather similar to a four-state Blume-Emery-Griffiths model. 

It is a well-known feature of quadrupolar ordering that its phenomenol- 
ogy closely imitates that of antiferromagnetic transitions, though the un- 
derlying order parameter is non-magnetic. The phase diagram in the H—T 
plane (H: magnetic field) was mapped in [11]. The salient features are the 
following: A sufficiently strong field applied in any direction will suppress 
AFQ ordering completely. On the phase boundary, a low-field regime of 
continuous transitions is separated by a tricritical point (H* « 2Tesla, 
T* ~ 4K) from the high-field regime of first-order transitions. This change 
in the character of the phase transition is shown in the field dependence 
of the specific heat. The nature of the magnetization curve changes dras- 
tically at T*. For T < T*, there is a steplike metamagnetic transition 
corresponding to the first-order transition from the low-T ordered phase to 
the disordered phase. For T* < T < T tr (H = 0) « 6.5K, there is a kind of 
a smooth metamagnetic transition, where the system crosses the second- 
order part of the phase boundary. For T > T tT (H = 0), the magnetization 
curve is completely smooth. We will show that a mean field treatment of 
the AFQ transition in the Ti-I^ scheme accounts for these observations 
quite well. 

In the absence of an external magnetic field, T3 and T4 type order 
parameters (i.e., quadrupolar moment and magnetization) are decoupled 
because the former is invariant under time reversal, while the latter changes 
sign. Switching on the magnetic field breaks time reversal invariance, allow- 
ing that quadrupolar moment and magnetization get mixed. We can also 
say that magnetic field, though it couples directly to the angular momentum 
J, may also induce quadrupolar moment. 

This is best illustrated by looking at the matrix which contains the 
matrix elements of the crystal field, the Zeeman energy for a field in the 
x-direction —h x J x , and also a term containing the quadrupolar moment 
\G\, within the basis (l)-(4) 



M(X,h x ,A) = 




(7) 



The field couples the singlet ground state to the O^-moment bearing 
excited state |r^ ). Therefore in the presence of a magnetic field, uniform 
quadrupolar moment is no longer "spontaneous". If at H = 0, we had 
to do with a transition to a ferro quadrupolar state, it would be smeared 
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Figure 1. The temperature dependence of the magnetic susceptibility. The onset of 
02-type quadrupolar order makes the two H \\ z and H \\ x behavior inequivalent. 

out in H ^ 0, and we no longer had a phase boundary to speak about 5 . 
However, for antiferroquadrupolar coupling, the appearance of the stag- 
gered quadrupolar moment is still symmetry breaking, and therefore a 
sharp phase transition remains possible also in an external magnetic field. 
Therefore, if we had no other evidence than that PrFe4Pi2 has sharp 
phase transitions in external magnetic field, and we adopted the ri-T4 
scheme, we would have to conclude that the ordered state could not be 
ferroquadrupolar, but only antiferroquadrupolar. 

2.3. MEAN FIELD RESULTS 

Doing the mean field theory [3] for AFQ ordering involves diagonalizing 
matrices like (7), and we do not give the details here. 

Fig. 1 illustrates that an AFQ transition, though of non-magnetic na- 
ture, may give a susceptibility which looks very much like what you expect 
from an antiferromagnet. The cubic (001) and (100) directions are equiva- 
lent in the para phase but the appearance of O^-type AFQ order makes the 
x-field susceptibility appear as "transverse", while the z-field susceptibility 
looks "longitudinal". Of course, instead of 0\ = J% — Jy we might have 
chosen Jy — J| or j| — J%, so in a crystal one would expect an equal 
mixture of the corresponding AFQ domains, and the susceptibility suitably 
averaged. The experiments may be taken to correspond to this. 

5 Purely in symmetry terms: applying a field in one of the cubic (100) directions, the 
symmetry would be lowered to Cah, and the decomposition of T4 of Oh in terms of the 
irreps of C^h would contain the identity which also comes from Ti ground state. 
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Figure 2. Left: The boundary of the antiferroquadrupolar phase in the H-T plane. 
The curve is drawn in black for first-order transitions, and in grey for continuous phase 
transitions. Arrow indicates the tricritical point. Right: The T-dependence of the tem- 
perature coefficient of the specific heat for H = 0, 2.5, and 4Tesla (in order of decreasing 
transition temperatures) . 



Fig. 2 (left) gives the phase diagram. The crystal field splitting A and 
the quadrupolar coupling A were chosen so as to get at least a rough numeri- 
cal agreement with the experimental phase diagram [11]. There is still some 
freedom in the parameters, but we found that a rather low A ~ 4— 6K has to 
be chosen (with z = 8, A 0.08/cb), if we want to get both T tT (H = 0) and 
the tricritical point right. These estimates are likely to be subject of some 
revision when further (especally dipolar) couplings are allowed for. - In 
spite of an overall similarity to the phase diagram based on experiments, we 
note that the low-T, high-i/ upcurving part of our present phase boundary 
represents a deviation, the reason for which remains to be clarified. 

The changeover to a regime of first-order transitions in higher fields is 
evident in the field dependence of the specific heat; the curves shown in 
Fig. 2 (right) bear a close resemblance to the measured ones. The same is 
true of the magnetization curves (Fig. 3) where we see a change from the 
regime of sharp metamagnetic transitions at low temperatures to continu- 
ous phase transitions at intermediate T's, and eventually smooth behavior 
in the para phase. 
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Figure 3. Antiferroquadrupolar order underlies the sharp metamagnetic transition at 
low T (left). With increasing T, we pass through second order transitions (middle) to 
smooth behavior (right). 

was supported by the grants OTKA T038162, and AKP 2000-123 2,2. 



References 

1. CD. Batista et al., Phys. Rev. B65, 180042(R), 2002, and references therein. 

2. A.T. Boothroyd, J. Alloys and Compounds 303-304, 489, 2000. 

3. A. Kiss and P. Fazekas, to be published. 

4. R. Shiina et al., J. Phys. Soc. Japan 66, 1741, 1997. 

5. A. Takahashi, H. Shiba, J. Phys. Soc. Japan 69, 3328, 2000. 

6. G. Mihaly et al., Phys. Rev. B61, R7381, 2000; G. Mihaly et al., at this conference. 

7. K. Penc et al., to be published. 

8. P. Fazekas, Found. Phys. 30, 1999, 2000. 

9. P.W. Anderson, Mat. Res. Bull. 8, 173, 1973; P. Fazekas and P.W. Anderson, Phil. 
Mag. 30, 423, 1974. 

10. K. Penc et al., at this conference. 

11. Y. Aoki et al., Phys. Rev. B 65 064446, 2002; J. Phys. Chem. Sol. 63, 1201, 2002. 

12. S.H. Curnoe et al, J. Phys. Chem. Sol. 63, 1207, 2002; Physica B 312-313, 837, 
2002. 

13. C. Sekine et al., Phys. Rev. Lett. 79, 3218, 1997. 

14. E.D. Bauer et al., Phys. Rev. B 65, 100506(R), 2002; K. Miyake et al., unpublished; 
D.E. MacLaughlin et al., Phys. Rev. Lett. 89, 157001, 2002. 

15. K. Iwasa et al., Physica B 312-313, 834, 2002. 

16. Y. Aoki et al., J. Phys. Soc. Japan 71, 2098, 2002. 



fazekas_hvar_b.tex; 1/02/2008; 23:09; p. 8 



